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ABSTRACT 

We present our optical observations of Swift GRB 070518 afterglow obtained at the 0.8- 
m Tsinghua University-National Astronomical Observatory of China telescope (TNT) 
at Xinglong Observatory. Our follow-up observations were performed from 512 sec 
after the burst trigger. With the upper limit of rcdshift~0.7, GRB 070518 is found 
to be an optically dim burst. The spectra indices (3 ox of optical to X-ray are slightly 
larger than 0.5, which implies the burst might be a dark burst. The extinction Ay 
of the host galaxy is 3.2 mag inferred from the X-ray hydrogen column density with 
Galactic extinction law, and 0.3 mag with SMC extinction law. Also, it is similar to 
three other low-rcdshift optically dim bursts, which belong to XRR or XRF, and mid- 
term duration(T 90 < 10, except for GRB 070419A, T 90 =116s). Moreover, its R band 
afterglow flux is well fitted by a single power-law with an index of 0.87. The optical 
afterglow and the X-ray afterglow in the normal segment might have the same mecha- 
nism, as they are consistent with the prediction of the classical external shock model. 
Besides, GRB 070518 agrees with Amati relation under reasonable assumptions. The 
Ghirlanda relation is also tested with the burst. 

Key words: gamma rays:bursts-gamma rays: observations-individual: GRB 070518. 



1 INTRODUCTION optical-to-X-ray spectral index /3 ox < 0.5), according to the 

„. . . . , fireball model, 
since the first optical counterpart of gamma-ray burst 

(GRB) was identified on 1997 February 28, more and more However, GRBs with an optical afterglow have been 
GRBs have been optically detected, especially after Swift studied widehr and in de P th ( e -S- GRB 060218, GRB 
was launched successfully, which allows early follow-up ob- 080319B). This has resulted in an increasing number of 
servations with ground-based optical telescopes. However, GRBs wlth known redshift, which has allowed statistical 
about 40 per cent of Swift GRBs still failed to be identified studies to be carried out ' The °P tical luminosity distribution 
with optical or infrared counterparts (Burrows et al. 2008). of GRBs vanes wldel y U P to several magnitudes (Kann et 
For the nature of the dark burst, several models have been aL 2006 ; 2008 )- After bein S transformed into a common red- 
proposed: extinction from dusty opaque star-forming regions smft ( e -S'> z = X )> the °P tlcal luminosity shows a bi-model 
( e.g. Groot et al. 1998; Reichart & Price 2002); high red- phenomenon (Nardini et al. 2006; Liang et al. 2006; Kann et 
shift, which diminishes the GRB optical luminosity (Lamb aL 2006 ; 2008 ) : °Pfacal]y luminous and optically dim. The 
& Reichart 2000), as the lya forest is shifted into the obser- bi-model distribution is also found to exist m other energy 
vational energy range; intrinsically dim bursts (Fynbo et al. bands (Gendre et al. 2008a,b). 

2001; Rol et al. 2005 ); the observational effect, As a result Liang & Zhang (2006) found that the apparent bimodal- 

of the late response to the bursts. Jakobsson et al. (2004) ity cannot be interpreted as a manifestation of the extinction 

have proposed an operational definition of dark burst (i.e. effect, and they suggested that there might be two types of 

progenitors or two types of explosion mechanisms in opera- 
tion. Considering the dark burst is usually a type of GRBs 

* email: xlp@bao.ac.cn with a faint optical afterglow, optically dim bursts might be 
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related to dark bursts (Nardini et al. 2006). However, the 
number of well-studied optically dim bursts is very limited. 
Thus, increasing the number of known optically dim bursts 
may help us to understand the origin of the bi-model distri- 
bution, the nature of optically dim bursts and dark bursts. 

One such interesting case is GRB 070518. It was de- 
tected by Gamma-ray Burst Alert Telescope (BAT) onboard 
Swift at 14:26:21 (UT) on 2007 May 18. The X-ray telescope 
(XRT) and ultraviolent/optical telescope (UVOT) on board 
Swift observed the counterpart beginning at 70 and 100 sec 
after the trigger, respectively. T90 (15-350 KeV) is 5.5±0.2 s 
,T5o is 2.9s (Sakamoto et al. 2008a). The values of duration 
mean the burst belongs to the short tail of the long group 
(Kouveliotou et al. 1993). Thus, it might be a classical long- 
duration GRB or an intermediate-long GRB (Horvath et al. 
2008). The ratio of fluence between the 25-50 KeV and 50- 
100 KeV band (Sakamoto et al. 2008a) is about 1.06, which 
makes GRB 070518 belong to an X-ray rich burst (XRR), 
according to the criterion|j given by Sakamoto et al. (2008b). 
The first given magnitude and coordinates of the counter- 
part are about 18 mag in white band reported by UVOT and 
RA(J2000) = 16h 56m 47.7s, Dec(J2000) = 55d 17m 42.3s 
(radius, 90 pen cent containment), respectively (Guidorzi et 
al. 2007) . The afterglow was also observed by other ground- 
based telescopes later. Besides, the redshift of GRB 070518 
was reported to be lower than 0.7 (Cucchiara et al. 2007) 
based on the photometry of the Swift UVOT. No other re- 
port about redshift of GRB 070518 has been presented in 
the literature up to now. In our analysis, the redshift of 0.7 
are used to investigate the properties of GRB 070518. 

In this paper, we report on optical photometric follow- 
up observations of GRB 070518 at Xinglong Observatory of 
National Astronomical Observatories, Chinese Academy of 
Sciences. Afterglow observations and results are described 
in §2. An anaysis and discussion are presented in §3, and a 
summary and conclusion are given in §4. In our calculations, 
a flat universe is assumed with matter density Qm = 0.3, 
cosmological constant = 0.7, and Hubble constant Ho = 
70 km s -1 Mpc -1 . The formula / orf~ c V~ /3 is used for the 
afterglow decay analysis. 



2 TNT OBSERVATIONS AND DATA 
REDUCTIONS 

We carried out a follow-up observation programme of Swift 
GRBs with the 0.8-m Tsinghua University-National Astro- 
nomical observatory of China optical telescope (TNT) at 
Xinglong Observatory, under the framework of the East- 
Asia Follow-up Observation Network (EAFON; Urata et al. 
2003;Urata et al. 2005). The TNT telescope is equipped with 
a PI 1300x1340 CCD and filters in the standard Johnson- 
Bessel system. The field of view of the TNT is 11.4x11.4 
arcmin, yielding 0.5 arcsec of the pixel scale. A custom- 
designed automation system has been developed for the 
GRB follow-up observations (Zheng et al. 2008). 



The follow-up observations of GRB 070518 were per- 
formed with the TNT from 512s after the initial Swift trig- 
ger, and lasted for about 5.5 hours until the dawn time. A 
sequence of white, R, V, and I bands images was obtained. 
After preliminary analysis, we first confirmed the counter- 
part of GRB 070518 (Xin et al. 2007) after the report of 
Swift UVOT. 

Data reduction was carried out following standard rou- 
tine in IRAF0 package, including bias and flat-field correc- 
tions. Dark correction was not performed, as the tempera- 
ture of our CCD was cooled down to —110 °C. Point spread 
function (PSF) photometry was applied via DAOPHOT task 
in IRAF package to obtain the instrumental magnitudes. 
During the reduction, some frames were combined in order 
to increase signal-to-noise (S/N) ratio. Flux calibration was 
performed by re-observing the field of GRB 070518, together 
with Landolt photometric standard stars with the TNT tele- 
scope on 2007 November 14 and December 12. Considering 
that the white band is not a standard band, we treated the 
lii/iite-band filter as a very broad filter. When we calibrated 
the white band data to R band magnitude, no significant 
difference was found in the colour terms, and the response 
functions of the two filters were very similar with about 0.07 
mag uncertainty. Since the redshift of GRB 070518 is no 
larger than 1, the white band data was calibrated to R band 
magnitude, for simplicity. The transformation uncertainty of 
0.07 mag is added to the final results. 

In addition, other optical photometries for GRB 070518 
are collected from the GRB Coordinate Network (GCN), in 
order to complete the optical light curves. All of these optical 
data are presented in Table 1. 



3 ANALYSIS AND DISCUSSION 

3.1 Light Curves and Chromatic Decay 

The Multiband light curves of GRB 070518 are shown in 
Fig. 1. Galactic extinction (Schlegel et al. 1998) is cor- 
rected for all V,R and I bands data, which corresponds to 
a reddening of Eb-v = 0.017 mag in the direction of the 
burst. The R- band data could be corrected for the flux con- 
tamination of its host galaxy, whose magnitude is estimated 
to be R = 23.8 ±0.1 mag (Garnavich et al. 2007b). The con- 
tamination is also interpreted as a combination of its host 
galaxy and supernova bump (Dai et al. 2008). As shown in 
Fig. 1, the R band light curve corrected for the flux con- 
tamination is labeled with "after", while the R band light 
curve without the correction is labeled with "before". The 
TNT data before 10 4 s are not affected by the contamina- 
tion very much. Moreover, the "after" R band light curve of 
whole data (including TNT data and GCN data) could be 
well fitted with a simple power law with a decay index of 
a a = 0.87 ± 0.02 (x 2 /dof = 54.6/17). The result is consis- 
tent with the fitting index of ~0.81±0.08 (xV do / = 14/10) 
only for TNT data. The general trends of light curves for 
the R (labeled "before") and / bands are similar; however, 



1 S(25-50KeV)/S(50-100KeV)=<0.72 C-GRB; 
0.72<S(25-50KeV)/S(50-100KoV)=<1.31 XRR; 
S(25-50KeV)/S(50-100KeV)>1.32 XRF. 



2 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under cooperative agreement with NSF. 

3 http:/ /irsa. ipac.caltech.edu/applications/DUST/ 
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Table 1. A log for the optical afterglow photometry of GRB 
070518. The reference time TO is 14:26:21 (UT) on 2007 May 18. 
The TNT white band (W) magnitudes were calibrated against 
the R band with the uncertainty of 0.07 magnitude. No data 
were corrected for the Galactic extinction of Eg^Y = 0.017. The 
flux contamination and extinction of host galaxy were also not 
corrected. 



T-T0 


Exposure 


Band 


Mag 


Telescope 


(mid day) 


(Sec) 






or Reference 




0.00636 


2x20 


W 


19.03±0.10 


TNT 


0.00699 


2x20 


W 


19.23±0.10 


TNT 


0.00804 


3x20 


W 


19.37±0.12 


TNT 


0.00883 


3x20 


W 


19.27±0.08 


TNT 


0.00961 


3x20 


W 


19.68±0.14 


TNT 


0.01053 


4x20 


W 


19.56±0.12 


TNT 


0.01157 


4x20 


W 


19.74±0.13 


TNT 


0.01275 


5x20 


W 


20.08±0.15 


TNT 


0.01406 


5x20 


W 


20.06±0.18 


TNT 


0.02446 


1x300 


R 


20.20±0.10 


TNT 


0.02814 


1x300 


R 


20.76±0.16 


TNT 


0.03723 


3x300 


R 


20.48±0.10 


TNT 


0.01237 




R 


19.5 


GCN6443 


0.23337 




R 


22 


GCN6483 


0.35948 




R 


22.2 


GCN6418 


0.37753 




Rc 


22.3±0.2 


GCN6458 


0.71837 




R 


23.03±0.05 


GCN6462 


1.45833 




R 


23.3 


GCN6426 


1.70837 




R 


23.36±0.05 


GCN6462 


3.60837 




R 


23.56±0.05 


GCN6462 


0.02065 


1x300 


V 


20.74±0.15 


TNT 


0.06383 


10x300 


1 


20.47±0.17 


TNT 


0.10197 


10x300 


1 


20.92±0.13 


TNT 


0.15271 


14x300 


1 


21.37±0.14 


TNT 


0.20660 


15x300 


1 


21.23±0.24 


TNT 




1000 1"* i° s 

Time after the burst (sec) 



to 8 to 7 



Figure 1. Multiband light curves of GRB 070518. Chromatic de- 
cay occurs between the X-ray and optical light curves. The energy 
range of the X-ray flux is 0.3 — 10 KcV. All R, V and / band data 
have been corrected for the Galactic extinction. Solid and open 
triangles present the R band data before and after the correction 
of the flux contamination of the host galaxy, respectively. The 
insert shows the flares in the X-ray light curve at the beginning. 



as a shallow decay. Therefore, the X-ray light curve displays 
a "steep-shallow-normal" decays (Zhang et al. 2006). 



the fluctuations exist in both light curves, which might be a 
result of a low S/N ratio. 

In order to investigate the decay properties between op- 
tical and X-ray light curves, we use the X-ray afterglow data 
downloaded from the website of the UK Swift science data 
center 0(labeled as USS) (Evans et al. 2007), where some 
analysis results are given, including the spectral indices for 
each segments in the X-ray light curve. These spectral in- 
dices are used directly for our later analysis. The X-ray light 
curve is also plotted by squares in Fig. 1, showing several 
drastic flares from the beginning up to about 250s after the 
burst trigger. Considering the sparse data between 2000 and 

3 x 10 4 sec, we fit the light curve after 250 sec by separat- 
ing it into two parts, one is between 250 sec and 2000 sec, 
and the other is from 3 x 10 4 sec to over 6 x 10 6 sec. The 
first part is well fitted by a single power law with an in- 
dex of a x ,i = 1.78 ± 0.11 (x 2 /dof = 0.76/3), while the 
second is well fitted by a single power law with a slope of 
a x ,i = 0.99 ± 0.08 (x 2 /dof = 5.83/8). As shown in Fig. 
1, the extrapolation of the first fitting is below the data 
detected after 1000 sec, suggesting that there is energy in- 
jection (Zhang et al. 2006) between the two segments, shown 

4 http://www.swift.ac.uk/xrt_curves/00279592/ 



3.1.1 First Decay in X-ray Light Curve 

The first steep decay in the X-ray afterglow light curve 
is usually explained by the Curvature Effect (e.g. Kumar 
& Panaitescu 2000; Zhang et al. 2009). According to their 
model, temporal index a and spectral index f) have a tight 
relation a = 2 + /3. For GRB 070518, the first slope a x ,i 
(1.78) of the X-ray light curve is definitely lower than the 
value (3.22, from 2 + /3) expected for high-latitude emission 
(Here, /3 = 1.22 from USS). However, the high-latitude ef- 
fect only provides an upper limit to the decay slope, as it 
assumes that the shell emission stops abruptly after the ini- 
tial pulse (Mangano et al. 2007). The decline could be made 
in two cases: one is that residual emission from the shocked 
shells is still present, in which, the high- latitude emission 
would contribute only a small fraction of the overall radia- 
tion; the other is that the X-ray band is below the cooling 
frequency, in which case, the expected decay slope of high- 
latitude radiation would naturally be shallower than 2. For 
GRB 070518, the first possibility might be more reasonable, 
from the analysis for other segments of the X-ray light curve 
in the later sections. 

3.1.2 Other Segments 

The X-ray afterglow is analyzed by applying the closure re- 
lations given by Zhang et al. (2006) to the X-ray spectral 
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and temporal indices. These relations help us to determine 
the location of the observed X-ray band relative to the syn- 
chrotron frequencies (v a ,Vm,Vc) and the environment where 
the burst occurred. 

The normal decay segment with a x ^ = 0.99 and f3 Xt 3 = 
1.11 ± 0.11, agrees with the non-injected closure relations: 
j3 = § and a = 3p 4 " 2 . These relations are consistent with two 
cases : when v x > v m in fast cooling phase or when v x > v c 
in slow cooling phase. Also, The environment condition [in- 
terstellar medium (ISM), Wind] could not be distinguished 
for the above cases. 

The application of the closure relations to the X-ray 
temporal and spectral index in normal decay segment, com- 
bining the analysis in § 3.1 for the X-ray afterglow during 
2000~2x 10 4 sec, suggests that there is energy injection be- 
fore the start of the normal decay phase. 

3.1.3 Flares in the X-ray Afterglow 

The X-ray light curve of GRB 070518 shows drastic flares 
between 80 and 200s since the trigger, as shown in the in- 
sert in the Fig. 1. Several GRBs were reported to have flares 
in X-ray afterglows (e.g. Bernardini et al., 2009; Krimm et 
al., 2007b). With the method of Krimm et al. (2007b), the 
ratios of raise time duration and the middle time for two 
big flares (peak times of about 106.9 s and 189 s, respec- 
tively) are calculated. AT/T for two flares was about 0.2 
and 0.1, respectively, and the AF/F was about 1. Such steep 
and large amplitude rises in the observed flux, as in GRB 
070518, could not be explained by a sudden increase in the 
external density (Nakar & Granot 2007). Another model, 
a patchy shell model, can produce the flare with AF/F of 
about 1 (Nakar & Oren 2004), athough it cannot interpret 
the flares with AT/T << 0.1. For two big flares of GRB 
070518, AT/T are about 0.2 and 0. ^respectively. Therefore, 
the model might be possible. Moreover, such dramatic flares 
are also related to refreshed shocks or the central engine ac- 
tivity (e.g. Zhang et al. 2006). This is more reasonable for 
the case of GRB 070518, if we consider the above analysis 
that energy injection was still occurred after these flares. 

3.1.4 Chromatic Decay and Afterglow Model 

Chromatic decay could be seen in the case of GRB 070518, as 
shown in Fig. 1. The X-ray afterglow shows a "steep-shallow- 
normal" decay (Zhang et al. 2005), while optical afterglow 
is decaying with a constant index of 0.87 (if the correct of 
the late contamination was right). The decay index of the 
optical light curve (a op t=0.87) is slightly smaller than that 
of the normal decay segment (0.99) of the X-ray afterglow. 
The normal decay segment is usually explained the external 
shock model (Sari et al. 1998), as expected for the optical 
afterglow. However, the origin of the normal phase might 
also be related to the internal energy dissipation (see as the 
discussion by Urata et al. 2007). Thus, the classical external 
shock model will tested with the burst. According to the 
model, the difference between the decay indices of the optical 
and X-ray afterglow should be between -0.25 for the uniform 
ISM case and 0.25 for the wind medium case (Urata et al. 
2007). For GRB 070518, a - a x = 0.87 - 0.99 = -0.12, 
ignoring the uncertainties of those values. The case of GRB 



070518 agrees with the prediction of the classical external 
shock model. Moreover, with the relation derived by Urata 
et al. (2007) : a„-a x = -0.25+s/(8— 2s), where s shows the 
condition of ambient matter density as a formate of n < r~ 3 , 
the value of s for GRB 070518 could be deduced to be 0.84. 
This suggests that the surrounding matter of the burst is a 
mix of ISM and wind-type medium. 

3.2 Testing the Amati and Ghirlanda relations 

The time-averaged spectrum of GRB 070518 from T-1.8s to 
T+4.5s is best fitted by a simple power-law model with a 
photon index F = 2.11 ±0.25 (Krimm et al. 2007a; Guidorzi 
et al. 2007), corresponding to the energy spectral index of 
1.11±0.25 {fi = r - 1). The fiuence in the 15-150 KeV band 
is 1.6±0.2xl0" 7 erg cm ~ 2 (Krimm et al. 2007a), which is 
lower than the average value of the Swift burst (Sakamoto 
et al. 2008a). 

3.2.1 Amati relation 

Fig. 2 shows the relation between E pea k,i in the rest frame 
and Ei SO (i.e. Amati relation; Amati et al. 2002). The data 
shown by open symbols and the best fitting parameters are 
obtained from Amati et al. (2008). For GRB 070518, spikes 
in the BAT light curve reported by Guidorzi et al. (2007) 
were only clearly seen under 50 KeV, which implies that 
E pea k,o might be less than 50 KeV. Besides, if E pea k.o is 
within the BAT energy range, the photon index T of a sim- 
ple power-law fit and E pea k, are well correlated with a re- 
lationship (Liang & Zhang 2005): 

logE peaki0 = (2.76 ± 0.07) - (3.61 ± 0.26)logF 
Sakamoto et al.(2009) also obtained a similar conclusion by 
combining the simulation: logE pea k = 3.258 — 0.829r 
With r = 2.11 and the two relations above, E pea k t0 = 
38.81!?; H an d 30.8 are obtained, respectively. The mean 
value of ~ 35 KeV is taken from the above calculations as 
its peak energy in the observational frame. Thus E pea k t i = 
(l + z) x E pea k t o ~ 60 KeV in the rest frame. As a result, 
GRB 070518 is consistent with Amati relation within 3 a 
uncertainties, as shown in Fig. 2. 

3.2.2 Ghirlanda relation 

Ghirlanda, Ghisellini & Lazzati (2004) proposed a relation 
between peak energy E pea k,i in the rest frame and the 
collimation-corrected energy E 1 : E pf , a k,i — 480 x (-E 7 ,5i) 0,7 
KeV, where -E 7 ,si = -E 7 /10 51 ergs and E y = (l~cos9)Ejj so . 
In this relation, the jet angle 9 is very important, which was 
given by Sari, Piran & Halpern (1999): 
9 = 0.161 x (t jet ,d/(l + z)) 3/8 x (n x 77 7 /£ 7 , ISO ,52) 1/8 , 
where z is the redshift, tj et ,d is the break time in days. For 
GRB 070518, the time of a possible jet break should be later 
than the last observation of X-ray afterglow, tj e t.d > 10 6 sec, 
according to the analysis of § 3.1. Thus, the opening angle of 
jet 9 should be larger than 25 degree, under an assumption 
of the typical values: n — 3 cm" 3 and r)- 1 — 0.2 (Ghirlanda et 
al. 2004). From the values of 9 and E 7 ,i SOj 52, -E 7 is estimated 
to be larger than 6xl0 49 ergs. 

However, E pea k t i has been estimated to be about 60 
KeV in the observation frame, as analyzed in §3.2.1. The 
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Figure 2. Correlation of E pea k % and Ei so ( the Amati relation). 
The solid triangle is shown for GRB 070518 when z = 0.7 and 
Ep ea fe,o = 35 KeV. Ei so is calculated using the Band function 
(Band ct al. 1993), assuming photon index Y2 = 2.11 of the simple 
power law is the higher photon index, and fixing Ti = 1 as the 
lower photon index. The data shown as open symbols are taken 
from Amati et al. (2008). 

E pea k,i gives an Ej of about 5.1xl0 49 , based on the 
Ghirlanda relation. This value is slightly less than that 
from the above analysis (6xl0 49 ergs), but the values are 
consistent with each other, considering all the assumptions 
we have made. Therefore, GRB 070518 might agree with 
Ghirlanda relation from above analysis. 



Figure 3. Comparison of R band luminosity of GRB 070518 with 
three other long bursts after transformation into the same red- 
shift of z = 1. The comparison long bursts are GRB 050416A, 
GRB 060512A and GRB 070419A, as labeled in the plot. GRB 
070518 was shown as a solid line based on the redshift of 2 = 0.7. 
The reference for data are as follows: GRB 050416A (Soderberg 
et al. 2007); GRB 060512A (Mundell ct al. GCN5118; Male- 
sani et al. GCN5122; Cenko GCN5125; Milne GCN5127); GRB 
070419A (Cenko et al. GCN6306; Williams et al. GCN6328; 
Milne et al. GCN6341; Nishiura ct al. GCN6308; Updike ct al. 
GCN6317; Lizuka ct al. GCN6316; Fu et al.GCN6311;Pozanenko 
et al. GCN6407; Garnavich et al.GCN6406; Hill et al. GCN6486). 



3.3 One of the Optically Dimmest Long Bursts to 
Date 

Nardini et al.(2006) found that GRB optical luminosity at 
12 hours (at rest-frame) after the trigger shows a bimodal 
distribution. With a redshift of 0.7 for GRB 070518, 12 h 
at the rest frame corresponds to 20.4 hours in the observer 
frame. Following the method of Nardini et al. (2006), and 
assuming /3 = 0.8, L opt ,i2h,i = 2 x 10 28 erg s _1 is obtained 
based on the redshift of 0.7, and a = 0.87. The low value 
of L op t : i2h,i °f GRB 070518 means that the burst is one of 
the optically dimmest GRB afterglows so far (Nardini et al. 
2008). 

We have also noted that the result might be affected 
by the subtraction of its contamination at late phase. In 
order to investigate this, we just calculate the optical lumi- 
nosity at the rest frame with the data without correction of 
the contamination of its possible host galaxy. In this case, 
L a pt,i2h,i ~ 6x 10 28 erg s _1 is obtained with the assumptions 
as above, indicating that the burst is really an optically dim 
burst. Moreover, if the real redshift of the burst was lower 
than 0.7, its luminosity would become much lower than the 
value given above. 



Another method to investigate the optically dim lumi- 
nosity of GRB 070518 is to compare its light curve with 
other long bursts directly. Kann et al. (2007) have collected 
GRBs in the pre- Swift and Swift era to compare the lumi- 
nosity of GRBs in the common redshift of 1. In this sample, 
GRB 050416A, GRB 060512 and GRB 070419A are at the 
lower edge of the distribution of optical light curves (Seen in 
the Figure 1 of Kann et al. 2008) . In order to make the com- 
parison, we have collected R band data of the above three 
dim bursts from GCN Circular and other literature. These 
data are transformed into a common redshift z = 1, with 
Popt = 1.3 for GRB 050416A (Soderberg et al. 2007) and 
fixed /3 pt = 0.8 for other bursts, GRB 070518, GRB 060512 
and GRB 070419A. In the transformation, all magnitudes 
are corrected for the Galactic extinction, and without con- 
sidering the extinction from their host galaxies. In addition, 
for GRB 070518, flux contribution of host galaxy is cor- 
rected. As shown in Fig. 3, the luminosity of GRB 070518 
is close to that of the three bursts, based on the redshift of 
0.7. Therefore, GRB 070518 should belong to optically dim 
bursts, like the other three bursts (e.g. GRB 050416A, GRB 
060512 and GRB 070419A). 
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3.4 Extinction 

One of the interpretations for the optically dim luminosity is 
the extinction and absorption from circumburst medium in 
its host galaxy. Here, we try to estimate the extinction along 
the line-of-sight to the burst. The observed neutral hydrogen 
column density along this line-of-sight to GRB 070518 de- 
termined by XRT onboard Swift is about Nh = (9.6 ±1.6) x 
10 20 cm~ 2 (Guidorzi et al. 2007), larger than the previously 
measured Galactic value (Nh = 2.2 x 10 20 ; Dickey & Lock- 
man, 1990). The excess value is about (7.4±1.6) x 10 20 cm~ 2 . 
If the excess value is from a previously unobserved over- 
density located within or close to our Galaxy, then the ex- 
cess extinction Ay = 0.64 would be obtained with the excess 
value of neutral hydrogen, based on the relation found by 
Predehl & Schmitt (1995): Ay = 0.56 x N H {W 21 cm,- 2 } + 0.23 
mag. If this is the case, extinction should not be a prominent 
cause for the optically dim luminosity. 

It is likely that the excess column density is located 
within the host galaxy of GRB 070518. As in the discussion 
for GRB 051022 (Nysewander et al. 2006), the rest frame 
Nh is given by the observed Nh,o x (1 + z) 2 ' 6 . The source 
frame is about N H ,i = (3 ± 0.6) x 10 21 with z=0.7. From 
the relation of host extinction Av,i = Av,obs x (1 + z) mag 
(Nysewander et al. 2006), Av,i = 3.2 mag would be given 
in the rest frame based on a redshift of 0.7. Clearly, under 
these conditions, high extinction would be a dominate cause 
of the optically dim luminosity of GRB 070518. 

However, all above estimates are made with an assump- 
tion of the Milky Way-type extinction. The optical extinc- 
tion from the measured Nh absorption should depend on the 
assumed dust-to-gas ratio. Some evidence shows that the ex- 
tinction law in the GRB host galaxy is likely to be a Small 
Magellanic Clouds (SMC) type. Using the small dust-to-gas 
ratio of the SMC (Pei 1992), the host extinction Av,i would 
be 0.3 mag when z = 0.7. The difference of the extinction 
value estimated above is very large, about a factor of 10. 
Actually, the real extinction law in the GRB host galaxy is 
uncertain in most cases, like GRB 050401, for example, the 
extinction Ay is estimated to be 0.62 ± 0.06 by fitting the 
optical spectrum, which is lower than the value Ay ~ 9.1 in- 
ferred from the X-ray column density (Watson et al. 2006). 
One reason for the case of GRB 050401 was that high energy 
release of gamma-ray and hard X-ray photons might destroy 
the circumburst dust in the GRB host galaxy (Watson et al. 
2006). The reason might also be the case of GRB 070518. 



3.5 Comparison with Other Three Optically Dim 
Bursts 

As shown in Fig.3, three Swift optically dim bursts GRB 
050416, GRB 060512A and GRB 070419A are labeled to 
compare with GRB 070518. In order to investigate the rea- 
son for the dim luminosity of GRB 070518, we try to com- 
pare them in different aspects. 

Table 2 presented the main parameters for the three 
bursts along with those of GRB 070518. There are several 
similar properties, as follows. 

(1) For all the bursts, Tgo is lower than 10s, except for GRB 
070419A. This means that they belong to the shortest part 
of the "long GRB" (Paciesas et al. 1999), or intermediate 
long GRB ( Horvath et al. 2008 ) classification. 



Table 2. A comparision of properties for four GRBs. Y is the 
photon index at 15-150 KeV. F is fluencee xlO -7 erg at 15-150 
keV. Z is redshift. 



GRB 


Tgo 


r 


F 


Z 


Ay 


XRR 


name 


00 










or XRF 


070518 


5.5 


2.11 


6.02 


0.7~1 




yes 


050416A' 


2.4 


3.4 


3.5 


0.6535 


0.2(0.32) 


yes 


060512 6 


8.6 


2.49 


2.3 


0.4428(2.1) 




yes 


070419A' 


116 


2.35 


5.6 


0.97 




yes 


a The 


references 


for 


GRB 


050416A arc 


Sakamoto 


ct 



al.(2005,2006), Holland ct al. (2007) and Kann et al. (2007). 
6 The references for GRB 060512 are Cummings et al. (2006), 
Starling R. et al. (2006), Bloom ct al. (2006) and Kann ct 
al.(2007). 

c The references for GRB 070419A are Cenko et al. (2007) and 
Stamatikos et al. (2007). 



(2) The photon indices of these bursts are larger than 2, 
which implies that the spectra are so soft that these bursts 
belong to XRRs and XRFs. 

(3) The redshifts of all these bursts are lower than 1. 

Another aspect concerns the absorption and extinction 
of the GRB host galaxy. Of the bursts mentioned above, only 
GRB 050416A has been studied extensively (e.g. Sakamoto 
et al. 2006; Mangano et al. 2007; Soderberg et al. 2007), 
and the extinction from its host galaxy has been estimated 
to be Ay = 0.2 (Holland et al. 2007), or A v = 0.32 ± 0.17 
(Kann et al. 2007). The two values are similar to the aver- 
age value (Ay=Q.2) of GRBs based on the study of Kann 
et al. (2007, 2008). However, they are significantly smaller 
than the extinction (Ay ~ 3.8) expected from the hydrogen 
column density inferred from X-ray observations of XRF 
050416A, assuming a dust-to-gas ratio similar to that found 
for the Milky Way (Holland et al. 2007). The low value of 
extinction from host galaxy of GRB 050416A implies that 
the extinction of host galaxy of GRB 070518 might be not 
large, only about Av = 0.3, as analyzed in §3.4. 



3.6 Optical to X-ray Energy Spectral Index 

Optically dim bursts might be related to dark bursts (Nar- 
dini et al. 2008). Considering that GRB 070518 is an opti- 
cally dim burst, we wonder whether the burst is a dark burst 
or not. 

The index of optical to X-ray spectrum /3 ox (< 0.5) at 
11 h after the burst trigger is used to define dark bursts 
(Jakobsson et al. 2004). Cenko et al. (2009) argued that the 
X-ray flux at about 11 hours might be contaminated by the 
shallow decay for some bursts. As a result ,the X-ray flux 
would be larger than the real X-ray flux radiated along with 
optical afterglow, which would make some normal bursts 
"dark" bursts. They calculated /3 ox for the bursts in their 
catalogue at 1000s after the trigger time. However, it is too 
early for most of the Swift bursts to carry out the calculation 
at 1000s after the trigger, because the X-ray light curves in 
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ing to the definition of Jakobsson et ai. (2004). The optical 
extinction in its host galaxy inferred from X-ray hydrogen 
column density is 3.2 with Galactic extinction law, and 0.3 
with SMC extinction law. 

Afterglow model has been applied to the light curves. 
Energy injection is found to occurr for the X-ray afterglow 
before the start of the normal decay segment. The behavior 
of the normal segment of the X-ray afterglow and the optical 
afterglow is consistent with the prediction of the classical 
external shock model (Urata et al. 2007), indicating that 
they come from the same origin. As an XRR burst, GRB 
070518 agrees with Amati relation. Moreover, it fills the gape 
of high-luminosity and low-luminosity bursts in the Amati 
relation, similar to GRB 050416A (Sakamoto et al. 2005). 
It is also found that GRB 070518 agrees with the Ghalanda 
relation within reasonable assumptions. 

Moreover, a comparison with three optically dim bursts 
(GRB 050416A, GRB 060512 and GRB 070419A and GRB 
070518) reveals that several similar properties were shared: 
burst duration T 90 <10s (expect for GRB 070419A, T 90 = 
116s), soft spectrum at BAT band T> 2, low redshift z <1, 
etc. 



Figure 4. The optical and X-ray flux density light curves. The 
optical corresponds to the R band, shown as open hexagonal. 
The X-ray data at 3 KeV is shown as solid squares. The spectral 
indices are calculated at different times (1000s, 11 hours and 10 5 s) 
by the interpolation of their light curves. The optical to X-ray 
spectral indices (3 ox are also shown. 

some Swift GRBs might be at the steep or shallow decay 
segment about this time. 

For GRB 070518, /3 ox at 1000s, 11 hours and 10 5 s after 
the burst are calculated by interpolation, as shown in Fig. 4. 
/3 ox for GRB 070518 is 0.666, 0.606 and 0.543 at the above 
three epochs, respectively. It seems that the burst is not a 
dark burst, for all of /3 ox are slightly larger than 0.5. How- 
ever, spectral index (5 x depends on the energy distribution 
p of electrons (Jakobsson et al. 2004), For GRB 070518, p 
should be 2.22±0.22, as j3 x = 1.11 ± 0.11 (from analysis of 
USS) and p — 2/3. In this case, the criterion for dark bursts 
should become /3 ox < 0.6 ± 0.1. Consequently, GRB 070518 
might be a dark burst (Jakobsson et al. 2004) or a gray burst 
(Zheng et al. 2009). 



4 SUMMARY AND CONCLUSION 

We have presented optical afterglow observations of GRB 
070518 from the TNT and an analysis of the optical and 
X-ray light curves. The optical afterglow shows a constant 
power law decay ~0.87, while the X-ray shows a "steep- 
shallow-normal" decay behavior. Based on the upper limit 
of redshift of 0.7, the optical luminosity at 12 hours in the 
rest frame was estimated to be 2 x 10 28 erg s _1 , which means 
that the burst is an optically dim burst. The conclusion is 
also supported by comparing with three other optically dim 
bursts: GRB 050416A, GRB 060512 and GRB 070419A. The 
spectral indices of /3 ox at 1000 s, 12 hours, 10 s s are calcu- 
lated. All /Sox were slightly larger than 0.5, indicating that 
GRB 070518 might be a dark burst or a gray burst, accord- 
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